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Abstract 

The  spheroidal  T-matrix  formalism  developed  by 
Hackman1-3  and  Sammelmann4'6  for  acoustic  scattering 
is  extended  to  electromagnetic  scattering  from  lossy 
dielectric  solids  in  a  conducting  medium.  The 
spheroidal  T-matrix  formalism  exhibits  superior 
performance  with  respect  to  the  spherical  T-matrix 
formalism  for  objects  that  deviate  appreciably  from  a 
spherical  shape.  Both  acoustic  (elastic)  and 
electromagnetic  scattering  are  solutions  of  the  vector 
Helmholtz  equation.  In  the  case  of  elastic  wave 
scattering,  the  displacement  field  has  3  degrees  of 
freedom  corresponding  to  the  2  polarization  states  of 
the  shear  wave  and  the  longitudinal  mode.  In  the  case 
of  electromagnetic  scattering,  the  electric  (magnetic) 
field  has  2  polarization  states  corresponding  to  left  and 
right-handed  photons,  but  lacks  a  longitudinal  mode. 

The  T-matrix  description  of  electromagnetism  mimics 
the  T-matrix  description  of  elastic  wave  scattering  in  the 
absence  of  a  longitudinal  mode.  Indeed,  the  stress 
tensor  of  the  displacement  is  replaced  by  the  exterior 
derivative  of  the  electric  field  in  Betti’s  identity  in  the 
derivation  of  the  T-matrix  formalism  of  scattering  from  a 
lossy  dielectric  solid.  Continuity  of  the  displacement 
and  surface  traction  is  replaced  by  continuity  of  the 
tangential  components  of  the  electric  and  magnetic 
fields  in  the  boundary  conditions.  In  the  case  of  a  time 
harmonic  field,  the  presence  of  a  finite  conductivity  in 
the  medium  is  represented  by  the  insertion  of  an 
imaginary  component  of  the  wavenumber  that  is 
proportional  to  the  conductivity  in  the  medium.  In  the 
case  of  complex  wavenumber,  the  Helmholtz  equation 
is  no  longer  a  self-adjoint  operator,  and  the  S-matrix  is 
no  longer  unitary.  This  article  describes  some  of  the 
features  unique  to  scattering  in  seawater  due  to  the 
large  conductivity  of  the  medium. 

1.0  INTRODUCTION 

This  article  describes  the  application  of  the 
spheroidal  T-matrix  formalism  developed  by  Hackman1'3 
and  Sammelmann4  6  to  electromagnetic  scattering  from 
isotropic  lossy  dielectric  solids  in  a  conducting  medium. 

Section  2.0  begins  with  an  introduction  to  the 
scalar  and  vector  spheroidal  wavefunctions,  followed  by 
a  derivation  of  the  generalization  of  Betti’s  third  identity 


for  electromagnetic  scattering,  and  concluding  with  a 
derivation  of  the  T-matrix  for  a  large  aspect  ratio 
dielectric  solid  and  a  perfectly  conducting  solid  based 
upon  the  generalization  of  Betti’s  identity.  Section  3 
concludes. 

2.0  SPHEROIDAL  T-MATRIX  FORMALISM 

The  spheroidal  T-matrix  description  of 
electromagnetic  scattering  from  a  conducting  lossy 
dielectric  solid  in  a  conducting  medium  is  derived  in  this 
section.  The  symbols  e  ,  y/+,  <j+,  and  s  ,  /./  ,  o  are 

used  to  denote  the  electric  susceptibility,  magnetic 
permeability,  and  conductivity  of  the  media  exterior  and 
interior  to  the  scatterer,  where  the  symbols 


k  =  J co2 s  li  +  icou  a 

+  v  ^  (21 

k  =  y  co  e  ji  +  ia)jU  cr 

are  the  wavenumber  in  the  exterior  and  interior  regions, 
respectively. 

The  notation  of  Morse  and  Feshbach7  is  used 
to  denote  the  regular  and  irregular  spheroidal  radial 
functions  of  the  third  kind  by  the  symbols  jeml (c  I  £)  and 

heml (c  I  £) ,  respectively.  These  functions  correspond  to 
the  radial  functions  of  the  first  and  third  kinds  denoted 
by  R“’  (c  1 and  R™(c  1 by  Flammer8  and 
Abramowitz  and  Stegun9.  The  symbols  Sml(c  1 77)  and 
S'2)(c  I  rj)  are  used  to  denote  the  angular  function  of  the 
first  and  second  kind.  Here,  c  =  kf  is  the 
dimensionless  wave  number,  where  /  is  the  semi-focal 

distance  of  the  spheroidal  coordinate  system.  The 
normalization  convention  of  Flammer8  for  the  expansion 

coefficients  d"1  (c)  of  the  angular  functions  of  the  first 
kind  is  adopted. 

The  regular  and  irregular  basis  functions  of  the 
scalar  Helmholtz  equation  in  the  exterior  and  interior  of 
the  scatterer  are  defined  as  follows. 
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The  angular  functions 

s2Sc±uh(P)  = 

I  s(m)  „  /  ,  [cos (m<p),cr  =  e 

J— , 
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are  the  spheroidal  equivalent  of  the  spherical 
harmonics.  Note  that  in  the  above  definition  of  the 
norm  of  the  angular  functions  the  integral  of  the  angular 
function  and  its  complex  conjugate  is  computed.  For 
complex  wavenumbers,  the  above  integral  without  the 
complex  conjugate  is  equal  to  a  unitary  matrix,  since 
the  angular  wave  equation  is  non-self-adjoint  in  the 
case  of  complex  wavenumbers,  and  the  corresponding 
eigenvalues  are  in  general  complex.  This  is  to  be 
contrasted  with  the  case  of  the  spherical  wave  functions 
where  the  angular  wave  equation  is  independent  of 
frequency  and  hence  self-adjoint  even  for  complex 
wavenumbers.  The  corresponding  radial  wave 
equation  is  superficially  non-self-adjoint,  but  it  is 
conformally  equivalent  to  a  wave  equation  that  is 
independent  of  frequency  and  self-adjoint  by  a  simple 
rescaling  of  the  radial  coordinate.  Hence  the 
eigenvalues  for  the  spherical  radial  wave  equation  are 
also  real.  This  difference  between  the  behavior  of  the 
spheroidal  and  spherical  wave  functions  for  complex 
frequencies  has  resulted  in  a  number  of  errors  in  the 
literature,  and  subtle  differences  in  various  identities 
and  expansions  for  real  and  complex  frequencies  in  the 
spheroidal  case.  For  example,  in  the  case  of  complex 
wavenumber,  the  scalar  Green’s  function  of  the 
Helmholtz  equations  is  no  longer  given  by  the  simple 
bilinear  summation  in  terms  of  the  regular  and  irregular 
spheroidal  wavefunctions. 

G(r,r')  =  expj+A  I  r  —  r  '  I]  / 4zr  I  r  —  r  '  I 
*  ik'Y'  w  .Ret//  , 

Rather,  it  is  given  by  a  bilinear  expansion  of  the 
following  form. 

G{r ,  r  ')  =  exp[+/£  I  r  —  r  '  I]  /  4zr  I  r  —  r  '  I 

=  ikT^,KR^r  (2'5 


The  matrix  C7".  is  a  symmetric,  unitary  matrix 

((£/  )”.  =  (E/*)”,),  which  reduces  to  the  identity  matrix 
in  the  limit  the  imaginary  component  of  the 
wavenumber,  vanishes,  that  is,  Im(c)  — >  0  . 

Define  the  following  basis  {V™,}  for  the  vector 

Helmholtz  equation  in  spheroidal  coordinates  in  the 
exterior  and  interior  regions  of  the  scatterer. 

C,=7-vx(0 

K 

k+ 

(2.6 

a  —  r 


f(/V,  0,+/: 2 77),  prolate 
(a,  a  ,a  )  =  *< 

’  l(/2!,0,-/277).  oblate 

The  first  pair  of  vector  basis  functions  is  a  pair  of 
solutions  of  the  transverse  vector  Helmholtz  equation, 
and  the  third  basis  function  is  a  solution  of  the 
longitudinal  vector  Helmholtz  equation.  Since  there  is 
no  longitudinal  component  of  the  electric  and  magnetic 
field,  the  vector  index  r  is  restricted  to  the  values  1 
and  2  corresponding  to  the  2  transverse  modes. 

Define  the  following  angular  functions,  as  the 
generalization  of  the  vector  spherical  harmonics. 

A<±),™,  =  W1-^2 dn  - 

V1  -  77 

Al±)2.™,  =  (n\F^d„+<P  r—dJsl I’  (2.7 

y/1 "  V 

A(±>  =  £s(±) 

3,  ami  J  ami 

The  above  functions  are  not  solutions  of  the  vector 
Helmholtz  equation.  They  satisfy  the  following 
orthogonality  condition,  where  the  upper  sign  is  used  for 
prolate  coordinates  and  the  lower  sign  for  oblate 
coordinates. 
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(2.8 


The  generalization  of  Betti’s  identity  for 
acoustics  scattering  to  electromagnetic  scattering 
follows. 

Theorem:  Let  u  and  v  be  a  pair  of  arbitrary  solutions  of 
the  source-free,  transverse,  vector  Helmholtz  equation 
in  the  region  V. 


-V  x  V  x  m  +  £2m  =  0  (2.11 

Define  the  exterior  derivative  T(u)t/  and  its  normal 
component  t(u)  as  follows. 

T(u)  .  =  V  u  —  V  u 

v  '  af)  a  p  pa 

(2.12 


t(u)  =naT(u) 

v  '  a  ^  '  era 


n<:'r  =  I  }  d  r/Sml  (c,  >i)Smr  ( c ,  r/) 

Aa i 


Here  h  is  the  outward  going  normal  to  the  boundary  of 
V.  Then  the  following  surface  integral  over  the 
boundary  of  V  vanishes. 

<^>dA{ti  •F(v)-F(m)»v}  (2.13 


The  matrix  Q(r  ”  is  a  symmetric  unitary  matrix  in  the 

case  of  complex  wavenumber.  In  the  case  of  real 
wavenumber  they  reduce  to  the  following  matrices. 

n!1."  =  a2-;  =  x  X 

l, l  l, l  ml  l 

2  +1 

C  r  2 

\lA,„A,,,r  (2-9 

Here,  Aml  are  the  eigenvalues  of  the  spheroidal  angular 
(radial)  wave  equation. 

At  this  point  it  is  worth  noting  that  the  vector 
basis  functions  {V™,}  have  the  following  asymptotic 
behavior  for  large  radial  coordinate  £  ->  °o  . 

Km  ->he*(c<&Ai.m 

(2.1 0 

C 

Km-t-defoJC'&Am 

c 

This  limit  coupled  with  Betti’s  identity  is  useful  in  proving 
various  orthogonality  relations  between  the  vector  basis 
functions. 


The  above  theorem  follows  from  Green’s 
theorem  and  the  fact  that  the  transverse  wave  equation 
can  be  recast  in  the  following  form  in  terms  of  the 
exterior  derivative. 

(— V  x  V  x  u  +  k'u)a  =  VaT(u)m  +  k~ua  =  0  (2.14 

Here,  the  exterior  derivative  T(u)  plays  the  role  of  the 

stress  tensor  in  elastic  wave  scattering,  and  the  interior 
derivative  t(u)  is  analogous  to  the  surface  traction  of 
the  stress  tensor. 

The  following  is  a  useful  corollary  of  the  above 
theorem. 

Corollary:  Let  S0  be  an  arbitrary  closed,  compact 
surface  about  the  origin,  the  sphere  at  infinity,  and  V 

the  volume  whose  boundary  consists  of  the  union  of 
these  two  surfaces.  Then  given  an  arbitrary  pair  of 
solutions  of  the  source-free,  transverse,  vector 
Helmholtz  equation  in  the  region  V,  the  following 

integrals  over  the  surfaces  S0  and  S X  are  equal. 

J dA{u  •  F(v)-F(m)»  v}  =|  dA{u  •  F(v) -  t(u)  •  v) 

s0  sx 

(2.15 

The  previous  theorem  and  its  corollary  allows 
one  to  write  the  following  orthogonality  relations  for  the 
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vector  basis  function  for  an  arbitrary  compact,  closed 
surface. 


(Nf)„  =  2X„,/n, 


JJ  clA{t{ Re  V  ^ )  •  Re  V.  ,  -  Re  V  ^  •  /(Re  ))  =  0 


=  +J  { Re  V <+)  •  t(E  )  -  f(Rey'+’)  •  E+ } 

5 


s 

\\dA{t{ ReV^J-V,.,..,,  -ReV^ 

S 

f-i / k(+i)‘~r >r  *  3 
[  0,  r  =  3  or  r  ’  =  3 


( Nra)n  = 

=  -}{y<+,./(£  )-/(y<+,).£  } 

5 


(V 


Tcrml,T'a'm'r 


=  —s;'s:s:'(+rrnTi: 

k 

+ 


(2.16 


jvr  =  +(V 


(2.18a 


(2.18b 


The  above  orthogonality  relationships  follow 
from  the  previous  corollary,  where  the  asymptotic  limit 
of  the  vector  basis  functions  {V  aml }  in  the  limit  |  ->  o°  , 
and  the  orthogonality  conditions  of  the  vector  spheroidal 
harmonics  {AiamI}  are  used  to  evaluate  the  surface 
integral  over  the  sphere  at  infinity. 

At  this  point  the  machinery  necessary  to  derive 
the  spheroidal  T-matrix  description  of  electromagnetic 
scattering  from  a  lossy  dielectric  solid  has  been 
introduced.  In  the  following,  the  indices  n  and  n’  are 
used  to  collectively  designate  the  indices  r,  ami  and 
z\a'm'V ,  respectively.  In  general,  the  index  r  is 
restricted  to  the  values  of  1  or  2,  since  the  exterior 
derivative  of  the  third  vector  basis  function  vanishes 
identically. 

The  incident  and  scattered  electric  fields  have 
the  following  expansions  in  terms  of  the  regular  and 
irregular  vector  basis  functions  in  the  region  exterior  to 
the  scatterer. 


E.nc  =  ^  a  Re  t// 1  ’ 


(2.17 


Here,  E  denotes  the  limit  as  the  exterior  electric  field 
approaches  the  surface,  and  E  denotes  the  limit  as  the 
interior  electric  field  approaches  the  surface. 

Next,  one  applies  Betti’s  identity  to  the  interior 
surface  fields.  Here  the  assumption  that  the  interior  is 
source  free  is  made,  and  the  interior  electric  field  can 
be  expanded  in  terms  of  the  regular  basis  functions.  In 
this  case  the  following  surface  integral  vanishes,  since 
Betti’s  identity  applied  to  a  pair  of  regular  vector  basis 
functions  vanishes. 

J{Re  <//;’•/(£)- /(Re  </’)•£}  =0  (2.19 

■S 


There  are  3  equations  for  the  interior  and 
exterior  surface  fields.  The  following  boundary 
conditions  are  utilized  to  rewrite  these  3  equations  in  a 
more  tractable  form. 


iix£  =  fix  E 

+ 

fix  H  =  hx  H 

+ 


(2.20 


Here,  the  source-free  Maxwell’s  equations  are  utilized 
to  relate  the  magnetic  field  to  the  curl  of  the  electric 
field. 


Using  the  above  orthogonality  conditions  for  the 
vector  basis  functions  and  the  above  expansions  for  the 
incident  and  scattered  fields,  one  arrives  at  the 
following  linear  equations  for  the  incident  and  scattered 
fields. 


1  -  - 

H= - Vx£  (2.21 

ia>ju 

Using  the  above  relationship  between  the 
electric  and  magnetic  field,  the  original  boundary 
conditions  can  be  rewritten  in  the  following  form. 
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h  x  £+  =  n  x  £ 

1  .  -  -  1  A  -  -  (2.22 
—  n  x  (V  x  £)  =  — n  x  (V  x  E  ) 

/£  +  £ 


Using  the  above  boundary  conditions,  coupled 
with  the  following  vector  identities, 


t  (£)  =  -n  x(Vx£) 

£  •  r(w)  =  — £  •  (n  x(Vxm)  =  -(£  x  n)  •  (  V  x  i?) 


(2.23 


one  can  re-express  equations  2.18  and  2.19  in  terms  of 
the  surface  fields  £  ,  and  t(EJ  . 

( Nf\  =  . 

=  +  J { Re y <+>  •t(E  )~  t( Re V(+> )  •  £  } 

5 


(N‘  a)  n  = 


(2.24 
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Re  Q,„ 

=  Jr(Rey,+,)»  A';’ 

5 

e„,„. 

• 

T 

V 

•4.., 

II 

ReM„, 

.  =  jRey<+)  .a;;’ 

s 

Mni 

II 

+ 

• 

>_ 

-  "+ 

£  = 

Jr(Rey,_))»  a';’ 

5 

P  ,= 

n,n 

1  n  n 

(2.27 


=  -J{V<+,.r(£  )-r(V<+,).£  } 

5 


The  solution  of  equations  2.26  for  the  scattered 
field  is  of  the  following  form. 


f=Ta  (2.28 

J  { Re  (//;’  •  f (£  )  —  -  r(Re  y/_) )•£_}  =  0 

s  T  is  the  spheroidal  T-matrix  given  below. 


The  surface  fields  £  ,  and  r(£  )  can  be  expressed  by 

the  following  expansion  in  terms  of  the  vector 
spheroidal  harmonics. 


£ 


I A 


(-) 


r(£  )  =  Xa,A!+) 


(2.25 


Substituting  these  expansions  into  equations 
2.24,  we  arrive  at  the  following  3  linear  algebraic 
equations. 


T  = -N' (ReQR'P -ReM)(QR‘P  -  M)' N  (2.29 

This  form  is  analogous  to  the  form  developed  by 
Hackman3  to  describe  the  acoustic  scattering  from  an 
elastic  solid. 

In  order  to  obtain  an  expression  for  the  interior 
field  in  terms  of  the  incident  field,  we  make  the  following 
expansion  of  the  interior  field  in  terms  of  the  regular 
basis  functions  in  the  interior  of  the  target. 

£,  =  ^ /;  Rc V,;  '  (2.30 


Nf  =  -  Re  Qa  +  ReM  J3 


Na  =  +Qa-  M  P 


(2.26 


0  =  Rcc-PP 


Here  the  matrices  are  defined  as  the  following 
surface  integrals,  where  there  are  no  longer  any 
unknowns  appearing  in  the  integrals. 


Next,  Betti’s  identity  is  used  to  obtain  the 
following  relationship  between  this  expansion  and  the 
surface  fields. 


5><-)„A.  =  J u(yl~') •  £  -vM  •  *(£_)} 

n'  s 


N 


(-) 


Tcrml  ,r '  a'  m’l’ 


k 


(2.31 
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Apply  the  boundary  condition, 

t(E)=—t(E+)  (2.32 

A+ 

to  equation  2.31  to  obtain  the  following  expression  for 
this  linear  equation  in  terms  of  the  surface  fields, 

E  and  t(E  ) . 

,  b.  =  VM  •t{E)}  (2.33 

;  a+  " 

Insert  the  expansions  given  in  equations  2.25 
for  the  surface  fields  into  equation  2.33  to  obtain  the 
following  linear  equation  for  the  interior  field. 

N^’b  -  Ba-C P 


Bnn,=  Jr(y<_)).  a';’  (2.34 


c  .  =  —  f  yl~’  •  a'V 

A+  s 

Solving  equations  2.26  and  2.34,  we  obtain  the 
following  expression  for  the  interior  field  in  terms  of  the 
incident  field. 

b  =  T~)a 


(2.35 

7,<")  =  (N{,y'(BR'P-C)(QR'P-M)'N 

In  the  case  of  a  perfectly  conducting  solid,  the 
boundary  condition  is  that  the  tangential  components  of 
the  electric  field  vanish  at  the  surface. 

h  x  E  =  fi  x  E  =0  (2.36 

The  discontinuity  in  the  tangential  component  of  the 
magnetic  field  is  proportional  to  the  induced  surface 
current.  In  this  case  the  T-matrix  for  a  perfectly 
conducting  solid  is  given  by  the  following  expression. 

T  =—N~l  ReM (M)  ' N  (2.37 

3.0  CONCLUSIONS 

The  spheroidal  T-matrix  for  a  lossy  dielectric 
solid  in  a  conducting  medium  was  derived  from  a 
generalization  of  Betti’s  identity  to  electromagnetic 
scattering.  The  resulting  T-matrix  is  formally  identical  to 


the  T-matrix  for  elastic  scattering  from  an  elastic  solid, 
where  the  stress  tensor  is  replaced  by  the  exterior 
derivative  of  the  electric  field. 

Future  work  will  focus  on  the  effects  of 
boundaries  on  the  scattering  from  proud,  buried,  and 
partially  buried  targets. 
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